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Abstract 

The electroelastic 4x4 Green’s function of a piezoelectric hexag¬ 
onal (transversely isotropic) infinitely extended medium is calculated 
explicitely in closed compact form (eqs. (73) ff. and (88) ff., respec¬ 
tively) by using residue calculation. The results can also be derived 
from Fredholm’s method [2]. In the case of vanishing piezoelectric 
coupling the derived Green’s function coincides with two well known 
results: Kroner’s expressions for the elastic Green’s function tensor [4] 
is reproduced and the electric part then coincides with the electric po¬ 
tential (solution of Poisson equation) which is caused by a unit point 
charge. 

The obtained electroelastic Green’s function is useful for the calcula¬ 
tion of the electroelastic Eshelby tensor [16]. 

Keywords. Electroelastic Green’s function, hexagonal medium, residues, 
piezoelectric materials 

Suggested PACS index: 77.22.-d 
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1 Introduction 


In the solution of many problems in physics Green’s functions or fundamen¬ 
tal solutions play an essential role. The elastic Green’s function tensor has 
been derived for isotropic media by Sir Thompson [1] and implicitely for 
anisotropic media by Fredholm [2], 

Especially the elastic Green’s function of the inhnite hexagonal medium has 
been calculated explicitely by Lifshitz and Rosenzweig [3] and by Kroner [4]. 

A Erst treatment of the electroelastic Green’s function has been given by 
Deeg [5]. Further contributions have been presented by Dunn [6], Wang 
[7] and Huang and Yu [8] in the framework of inclusion problems. But no 
compact explicit form for the electroelastic Green’s function has been given 
there. Thus, according to the author’s knowledge, a compact closed form 
representation of the electroelastic Green’s function does not exist. But, due 
to the widespread interest in piezoelectric materials a compact explicit rep¬ 
resentation is highly desirable. 

The residue method which is applied here is a useful tool to obtain Green’s 
functions in media with hexagonal symmetry as shown by Michelitsch [11] 
and Michelitsch and Wunderlin [12] for a treatment of the incompatibility 
problem. 

In the following we calculate using this residue method the electroelastic 4x4 
Green’s function for the hexagonal infinite piezoelectric medium. We shall 
obtain this Green’s function analytically in closed form (eqs. (73) ff. and 
(88) ff., respectively). To obtain a compact formulation the use of a con¬ 
venient symmetric tensor basis is essential [9, 10]. The basis tensors used 
here are very useful for the solution of several problems in hexagonal media 
[11, 13, 14]. 

The result for the electroelastic Green’s function presented here may be useful 
for the treatment of many problems, e.g. the inclusion problem in piezoelec¬ 
tric hexagonal (transversely isotropic) media. 

2 Basic Equations 

We start from the held equations for the stress tensor cr and the dielectric 
displacements D. The equilibrium conditions for the stresses are: 
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^j^ij •) ^ij ^ji (^) 

dj indicate the spacial derivatives and K the density of body forces, re¬ 
spectively. A further held equation describes the conservation of free electric 
charges: 


diDi = pe ( 2 ) 

Pe represents the density of free electric charges. The constitutive equa¬ 
tions (material law) of the piezoelectric medium which connect the electric 
held E and the elastic deformation e with the dielectric displacement D and 
the stress a are given by: 


Cijkl^^kl (^kijEk 

Di Gikl^^kl T ^ikEk 

Cijki, ^kij and pij denote the elastic moduli, the piezoelectric moduli 
and the dielectric moduli, respectively. They have the symmetry proper¬ 
ties C'ijki ^kiij ^jiki ^ijiki ^iki ^iik aud pij Vji' Iritroduciug the 
electric potential <h and the elastic displacement held u the ansatz for the 
strain e and the electric held E is given by: 

~ 2 + djUi) 

Putting ansatz (4) using (3) into the held equations (1) and (2) we obtain 
a 4 X 4 diherential equation of degree two for the held U = {u, <h) of the form: 


r{V)u + j^ = o (5) 

V indicates the gradient operator. Here we have introduced the gen¬ 
eralized force density ^ = {K,—pe). The symmetric 4x4 second order 
diherential operator (V) can be written in the form: 


r(v) 


T(V) t(V) 
t^(V) r(V) 


( 6 ) 


Here T (V) is a 3 x 3 tensor operator and represents the elastic part and 
is given by: 
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Tij (V) = Cipjqdpdq 


( 7 ) 


t (V) is a (3 X 1 tensor) vector operator given by 


ti(y) = Cpiqdpdq ( 8 ) 

and represents the piezoelectric coupling. 

Finally the (1x1 tensor) scalar operator 

r (V) = -Vpqdpdq (9) 

describes the dielectric part. 

The vector held lA can then be represented by the 4x4 electroelastic Green’s 
function ^ according to 

tl{r) = J Q (r — r'^ T (r'^ d^r' (10) 

r denotes the space point. The electroelastic Green’s function Q then is 
dehned according to 


r{V)gir) + 5^{r)l = 0 ( 11 ) 

^3 represents the three-dimensional (5-function and 1 denotes the 4x4 
unit matrix. 

The electroelastic Green’s function Qpq {p,q = 1,2, 3,4) has the following 
physical interpretation [6]: 

Qmj (r) (m,j = 1,2,3) is the elastic displacement at spacepoint r in the m- 
direction caused by a unit point force at spacepoint r = 0 in the j-direction; 
^m 4 (^) (^^ = 1, 2, 3) is the elastic displacement at spacepoint r in m-direction 
caused by a unit point charge at spacepoint r =0; 

(r) (j = 1, 2, 3) is the electric potential at spacepoint r caused by a unit 
point force at spacepoint r = 0 in the j-direction; 

^44 (r) is the electric potential at spacepoint r caused by a unit point charge 
at spacepoint r =0. 

We note that our electroelastic Green’s function dehned in equation (11) 
is symmetric, i. e. Opq = Qqp. This property follows from the symmetry of 
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the operator 7^^ (V) = Tsr (V) from equation ( 6 ). For our further calcula¬ 
tion we use the following convenient representation for the Green’s function 
[5, 9, 10, 11, 12]: 

( 12 ) 

This relation can be derived from (11) in straight forward manner by 
using Fourier transformation [11, 12], One gets T (^) from equation ( 6 ) by 
replacing di by {i = 1,2,3) in equations (7), ( 8 ) and (9). The vector ^ is 
then given by [ 11 , 12 ] 

^ (a) = ei cos a-f 62 sin a (13) 

The vectors form a useful orthonormal basis 


ei 


^ -y ^ 


/ —zx \ 


f X \ 


_ j_ 


^ _ 1 


X 

1 ^2 prf^ 

-zy 

1 ^3 f 

y 

V 0 ) 


\ P" ) 


{ ^ ) 


(14) 


where and r = re^. The orientation of this 

coordinate system can be expressed by e* = x e^. €ijk denotes the 

antisymmetric permutation tensor. 


3 Residue Calculation 

Our goal is to formulate the residue calculation ansatz to obtain the Green’s 
function from equation ( 12 ). 

First of all we introduce for convenience a useful orthonormal basis to rep¬ 
resent the 4x4 matrix 7^ (^) and later 7 (^): 



( \ 


( “^2 \ 


( 0 \ 

/ 

1 

6 

1 



0 


6 

0 

. 5 

0 

•j 

1 

, 64 = 


1 0 ) 


1 0 ) 


10 j 

V 


(15) 
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(■Cb = + ^ 2 ) = ■Cs)- This basis represents the hexagonal symmetry. 

Cfe and eb± are parallel to the basal plane and Cc represents the c-direction. 
64 comes into play because of the electric potential component <h. We then 
obtain the following useful representation; 


r(i) Th_i_6;5_i_ G^_i_ -j- TiqG}) Gb ^ ^ T T^Gq Gq 

+tb4 (Sf, ;8) 64 + 64 (g) Cfe) + tc4 (Gc 0 64 + 64 ® Cc) + TG4 0 64 

(16) 

The scalar quantities Tb_L, Tb, Tbc, T^, tb, 4 ,tc 4 , and r correspond to the 
tensors T, t and r from equations (7), ( 8 ) and (9), respectively (0 indicates 
dyadic multiplication). They are obtained as: 


Tb± = Cqq^I + €44^^, 

(17) 

Tb = Cii^^ + €44^^, 

(18) 

Tbc = (C13 + C44) ^h^c, 

(19) 

Tc = €44^1 + 

(20) 

th4 = (631 + 615) 

(21) 

tcA = 6l5^b + 633^0) 

(22) 

r = - (hii^b + V33C) 

(23) 


Cab — {Cii, C44,Cee,C13,C33} denote the elastic, Gia — {615,631,633} the 
piezoelectric, and rjij = {rjujVss} fhe dielectric moduli of the hexagonal 
material. Subscripts A, B represent Voigt’s notation whereas represent 
cartesian subscripts, respectively, 
r-^ (0 can be written as: 

T-' (i) = ^ (24) 
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Here Ajj denotes the matrix of 3 x 3 subdeterminants (multiplied by the 
prefactor (—and / (^) the determinant of 7~, respectively. 

We now can write for A by using (16): 

A (^) = Ab±eb± <S) eb± + A^Sb «) + Abe (e?, «) Cc + Cc ® Cb) + A^Cc ® Sc 

+Ab 4 (efo 0 64 + 64 (g) 6fe) + Ac 4 (6c 0 64 + 64 0 6c) + A464 0 64 

(25) 


Here the scalar quantities are introduced: 

Ab± = T (^TbTc — T^c) — (tc4^b — 2fc4tb4^bc + ^M^c) , (26) 

Ab = Tbx (TcT - 4) , (27) 

Abe = —T\,± (TbeT — tb4A4) , (28) 

Ae = Tb_L (TbT — tb4) ; (29) 

Ab4 = ^b± (^beA4 “ , (30) 

Ae4 = —^b± (^bte4 “ ^betb4) , (31) 

A 4 = Tb_L ^TbTe — T^c) (32) 

The determinant f {$,) oi'T then takes the form 

/(^) = Detr(^)=T(axAba(0 (33) 


Putting equations (17)-(23) into (33) shows that only terms proportional 
'Cb^'Cc”^” (n = 0,1, 2, 3,4) appear. This is a unique property of the hexagonal 
medium. As we shall show this is the central property which is needed to 
solve our integration problem (12) explicitely. Thus (26) is a polynomial of 
degree 3 in a when we put a = ^b/'Ce- We then can write 

Ab± = P(a)e® 
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where P (a) is a polynomial of degree 3 in a and takes the form: 


P (a) = Aa? + Ba? + Ca + D ( 35 ) 

We obtain for the coefficients A, B,C, D: 

A = —? 7 iiCiiC 44 — Ciie^5 ( 36 ) 

B = —^ 7331 ^ 11^44 — ■'711 (l^lll ^33 — 2C13C44 — C43) — C44ef^ — 2^11615633 

+2 (Ci3 + C44) 615 (631 + 615) — C44 (631 + 615) 

( 37 ) 

C = —7733 (C11C33 — 2C13C44 — C33) — 7711C33C44 — 2615633C44 — 633C11 


+2633 (631 + 615) (Ci3 + C44) — C33 (631 + 615)^ 

(38) 

D = — 7733 I 733 C 44 — 633 C 44 (39) 

Thus we can factorize the determinant / according to: 

/ (^) = + ^ 1 ) {a + A 2 ) {a + A 3 ) (a + A 4 ), (40) 

with T (a)^^_L = Cq 6 {a + Ai) (Ai = C 44 /CQQ) and 

P (a) = A (a + A 2 ) (fl + A 3 ) (a + A 4 ) (41) 

The A's from equations (26)-(32) yield with a = ^b/'Cc = 1: 

Ab± (a) = P (a) = Aa^ + Ba^ + Ca + D (42) 

The numbers A 2 , A 3 , A 4 are the zeros of the equation: 

Aa^ — Ba^ + Ca — D = 0 (43) 


with above coefficients A,B,C,D from equations (36)-(39). The Ai are ma¬ 
terial quantities and fully determined by the moduli C, e, 77. Furthermore, 
the subdeterminants (27)-(32) yield: 
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Ab (a) — — (Cgefl + C44) (rjiia + 7733) ((^440 + C33) + (e^^a + 633)^ , ( 44 ) 

Abe (fl) = \/fl (^^66® + C44) [(631 + 645) (CisO + 633) + (?7iia + 7733) (Ci3 + C44)] , 

( 45 ) 

Ac (a) = — (Ceea + C44) (7711a + 7733) (C^a + C44) + a (631 + 645)^ , ( 46 ) 

Ab4 (a) = \/a (Cgga + C44) [(C13 + C44) (e^^a + 633) — (C44a + C33) (634 + 645)], 

( 47 ) 

Ac 4 (a) = — (Ceea + C44) [{Cua + C44) (eisa + 633) — a (C43 + C44) (634 + eis)], 

( 48 ) 

A4 (a) = (Cge® + C44) a^C4iC44 + a {CiiC^‘4 — 2C13C44 — + C33C44 ( 49 ) 

The dependence on a is obtained by putting = i/a and .^c = 1 in equa¬ 

tions ( 17 )-( 23 ) and ( 26 )-( 32 ), respectively. 

To evaluate ( 12 ) we make use of the properties: 

A,, (|A) = A^^A,, (^) ( 50 ) 

A denotes an arbitrary scalar number. The components Aij are homoge¬ 
neous functions of degree 6. Whereas the determinant / is homogeneous of 
degree 8: 

/(|A) = AV(^) ( 51 ) 

Thus 

T-' («) = ( 52 ) 


10 



We now introduce the complex vector [12] 

7 (a) = (a) (53) 

with ^ (a) from equation (13) {i denotes the imaginary unit). (53) can then 
be written as 


7 (s) = h*s + h 

Here we have introduced the complex variable 

s = 

and the vector 


( 54 ) 

( 55 ) 


h = ei + ie2 (56) 

with the basis vectors ei 2 from equations (14). When putting A = 2e*" 
equation (52) then can be written as: 


r-^ (0 

Furthermore we observe 


4 

fhis)) 


(57) 


ds = 2isda (58) 

Then the integral (12) can be transformed into a complex integral over 
the unit circle: 


= 4A7^d£ 

SttV J\s\=i f (7 (s)) i 

Using the residue theorem we can rewrite integral (59) as 


(59) 


Qir) 


1 

STT^r 


27rf ^ Res 


( 4A (7(g)) \ 

w/ (7 (^)) ) 


(60) 


To evaluate (60) we have to hnd all zeros of / (7 (s)) which are located 
within the unit circle. To hnd these zeros we observe that /( 7 (s)) is a 
polynomial of degree 8 in s. Thus there exist 8 zeros Sj of /. We assume 
here that / has no multiple zeros Sj. 
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Following [12] we can conclude that there are 4 pairs of zeros si,si (/ = 
1, 2, 3,4) having the property: 


\sisi\ = 1 (61) 

There are only four zeros s/ which are located within the unit circle. Four 
zeros S] lie outside the unit circle. Thus the zeros can be written in the form 
[ 12 ]: 


Si = e2'^*e+2^h 


(62) 

(63) 


where 'ipi > 0 {I = 1,2, 3,4). Here the degenerate case |si| = |sj| is excluded 
by our assumption of no multible zeros. Only the residues corresponding 
to the four zeros (62) lying within the unit circle contribute to (60). The 
residues at s = s; are obtained as: 


Res 


( 4A (7 (g)) \ 

V*/ (t (^)) / 


Finally we arrive at 


4A (7 {si)) 

,d/ (7 (s))| 

1 Ls= 


( 64 ) 




47 rr df (7 (g)) , 
ds 

The zeros si of / (7 (s)) within the unit circle yield 


( 65 ) 


y/ Aip^ + z‘^ — r 
VAip^ + 2 :^ + r 

To obtain ( 66 ) equation (40) together with (54) and (56) is used. 


( 66 ) 


4 Explicit form of the Green’s function 


Evaluating (65) the Green’s function takes the form: 


g{r) 


4 

6 ^ VAip^ + ~ ^l) 


07 


( 67 ) 
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(/ = 1, 2, 3,4) are the solutions of the equations 


/(O = 0 ( 68 ) 

or 

■Cl + '^2 + = 0 (69) 

and 

= iix + i2y + 62 : = 0 (70) 

with ,^36 = 1 . Thns eqnations ( 68 ) and (70) have the solntions (/ = 1, 2, 3,4): 

-zx + iy^jAip^ + 2;2 (71) 

^ 2 ^ = -zy - ix^JAip^ + z'^ (72) 


{p = y/x^ + y'^). This formulation is in complete analogy to the resnlt 
obtained by Kroner for the Green’s fnnction tensor of the hexagonal elastic 
medinm [4] by using Fredholm’s method [2], 

The Green’s fnnction (67) then assnmes the form: 

Q (r) = (p, z) 60 ® + Qpp (p, z) + Qp^ (p, z) (e^ ®ez + e^® Sp) 

+Qzz (p, z) 6^0 6^ + Qpi (p, z) {Bp 0 64 + 64 0 Bp) 

+ QzA (P, z) {Bz 0 64 + 64 0 Bz) + ^44 (p, z) 64 0 64 

(73) 

Here we have introduced the following basis: 

( x\ /-p\ /0\ /0\ 

_ \ y _ \ X _0 _0 

~ P 0 ’’ ~ J ) g ■, ^z — , 64 — Q 

Vo/ VO/ Vo/ Vi/ 

(74) 

To obtain (73) we have used (25) together with (71) and (72). Eqnation 
(67) then yields the Q's introduced in (73) as: 
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1 ^ p^K{a =-Ai) ^ z‘^Ty,{a =-Ai) 

[p, Ar ^ 


47rACee ^ p^Aip'^ + z'^ n|=i,0¥0 


^pp (P) ^) 


E 


p^Ab± (g = -Ai) - z'^T^, (g = -Ai) 

AirACee^i p'^VAp^ + z'^ UUhim 

(7 ^ _ (~^) ^bc (q = ~A;) _ 

’ 47rAC66 pVA«p^ + 2 :^ nl=i,ovn (A 7 - Ai) 


GpA (p, ^) = 
Gzz (p, = 

^2:4 (P) 2 ^) 

^44 (P, 2 :) = 


E 


{-z) rb4 (a = -Ai) 


471 ACm pVAp^ + z'^ UUi,{j^i) - ^i) 


■E 


Ac (a = -Ai) 


471 ACm VAp’^ + z'^ 117=1,(j^z) (A - A) 

1 ^_ Ac4 (a = -Ai) _ 

47rAC66 VAATAuUpim (A - A) 


E 


A4 (a = -Ai) 


47rAC66 ^ VAp^ + n|=i,oyz) (Aj - A«) 
Here we have introduced the quantities 


(75) 

(76) 

(77) 

(78) 

(79) 

(80) 
(81) 


Abi (a) - Ab(a) = aFi, (o) (82) 

to arrive at 


rb(a) = (Cii - Cee) (Tc (a)r (a) - tc4 (a)) 

~ (^^13 + G 44 ) T (a) + 2 (C 13 + C 44 ) (csi + 645 ) tc4 A) (83) 

+ (631 + aA (fl) 


Abe (a) = A/ahbc (a) ( 84 ) 

where (compare eq. ( 45 )) 
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Tbc (fl) — (Ceefl + C 44 ) [(csi + 645 ) (ei^a + 633 ) + {r]iia + 7733 ) (C 13 + C 44 )], 

(85) 

and 


Ab4 (ci) = y/aTbila) (86) 

where (compare eq. (47)) 


rb 4 (a) — (CqqQ + C 44 ) [(Ci 3 + C 44 ) (CisO + 633 ) — (^440 + C 33 ) (631 + 645 )] 

(87) 

Equation (73) together with (75)-(81) represents the electroelastic Green’s 
function explicitely in compact form. Especially for the calculation of the 
electroelastic Eshelby tensor the use of representation (73) may be conve¬ 
nient. 


The cartesian representation of the the electroelastic Green’s function 
becomes: 


gw-E h , x 


^ 9il 9 i 2 9il 9u ^ 

9^1 922 £/S 92 I 

(i) (/) (/) (/) 

9i3 923 933 934 

I (0 (0 (0 (0 I 

V 9i4 924 934 944 ) 


Here we have used the abbreviations: 


(0 

^11 = TT 


P 


(0 _ 
922 — 




y'^z^ — {Aip^ + z'^) 


+ Ab± i—Ai) 


(0 ^ 
9i2 = TT 


(0 

^13 = 


-rb {-Ai, 

1 


xy {Aip^ + 2z^) 




-Ebc (-71/) 


P" 

xz 

p2 


( 88 ) 


(89) 

(90) 

(91) 

(92) 
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(0 

^23 

1 

~Il 

-hbci 

-Ai) 

1 1 

(93) 


(0 

^33 = 

1 , 

^ vrAci 
oi 

-Ai] 

1 

(94) 

(0 

£/14 

1 

~'Si 

-rb4( 

-Ai) 

xz 

p2_ 

(95) 

(0 

f/24 

1 

—rb4 ( 

-Ai) 

to 

1 _1 

(96) 


(0 

g\l = 

1 , 

^ Ac4 
01 

{-Ai 

) 

(97) 


(0 

gli = 

■ ^A4l 
C-l 

-Ai] 

) 

(98) 


and 

4 

Si = dTrCee^ {Aj — Ai) (99) 

i=i,OVO 

As we shall show in the appendix the Green’s function (88) yields in the 
case of vanishing piezoelectric coupling = 0 the well known results: For 
the elastic part Qij (i,j=l,2,3) Kroner’s elastic Green’s tensor of the hexagonal 
medium [4] and the dielectric part ^44 then represents the solution of Poisson 
equation of a unit point charge in a hexagonal dielectric medium, whereas 
the (j = 1, 2, 3) are vanishing. 

5 Appendix 

Here we consider the case of vanishing piezoelectric coupling {eijk = 0). As 
a consequence the terms (eqs. (21), (22)) 


tb4 = G4 = 0 (100) 

are vanishing. Thus the determinant / (eq. (33)) simplihes according to 
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( 101 ) 


/ (a) ^ DetT (a) ^ r (a) Tbi (a) (t^, (a) (a) - (a)) 

and can be written as 

/ (a) = —TjiiCiiCi 4 CQQ (a + Oi) (a + 02 ) (a + 03 ) (a + a 4 ) ( 102 ) 

where ai = Ai denote the zeros of / (—a) in the deconpled case. 

C44 

ai = -pr- 
‘-'66 

represents the zero of Tb_L(—a). The nnmbers a 2,3 are the zeros of the 
qnadratic eqnation 



Tb (—a) Tc (—a)—(~®) — CiiC 44 a^+(^C 43 + 2 C 13 C 44 — CiiC 33 ^ a+C 33 C 44 — 0 

(104) 

The zeros ai for / = 1, 2, 3 (elastic part) are those introdnced in [4], a 4 is 
given by 



(105) 


and represents the zero of r (—a). Furthermore we hnd when using (42)-(49) 
together with (82)-(87) 


Ab± (a — —a;) — —rjiiCiiCii (a 2 — a/) (a 3 — ai) — ai) (106) 


Ab (a — —a;) — —r/iiCee (C 33 — aiC^i) (ai — a;) (a 4 — a;) (107) 


Tb (a — —ai) — r]ii (a 4 — ai) (Cee — Cn) (C 33 — aiCii) + (C 13 + C 44 )^ (108) 


Tbc (a = -ai) = rjiiCm (C 13 + C44) (ai - a;) (a4 - ai) (109) 

(a = -ai) = -r/iiCee (C 44 - a/Cn) (ai - a;) (a 4 - a;) ( 110 ) 
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and 


T^i^a =-ai) = Aci{a =-ai) = 0 ( 111 ) 

to obtain Qj 4 = 0 (j = 1, 2, 3) and 

A4 (a = —ai) = C11C44C66 (cii — ai) (02 — ai) (03 — ai) (112) 

and for / = 1, 2, 3: 

Si =-r]ii{a4-ai)Ei (113) 

with 

3 

E; = 47rCiiC44C66 (dj — di) (114) 

i=i.0¥0 

The terms (114) coincide with Kroner’s (corrected) ”E;” (the terms E^ 
dehned in [4] have to be corrected by a prefactor —C11C44C66 to obtain the 
correct result for the elastic Green’s tensor [15]). 

For / = 4 we obtain; 


S4 = -Anr]iiCiiC44CQQ (oi - 04) (02 - 04) (03 - 04) = -47r?7iiA4 (-04) ( 115 ) 


We observe from equation (106)-(110) the properties: 


Ab± {d = —04) Ab (a = —04) Tb (a = —04) Tbc (n = —04) Ac (a = —04) ^ 

S4 S4 S4 S4 S4 

(116) 

(T4 7^ 0, eq. (115)). Thus the term / = 4 in sum (88) does not contribute 
to the elastic components Qij {i,j = 1, 2, 3). Using equations (112) and (113) 
together with (114) we obtain for / = 1, 2, 3: 


A4 (a = —di) 
Si 


0 


(117) 


Thus there are no contributions to sum (88) for / = 1, 2, 3 to the dielectric 
part ^44. For / = 4 we hnd from equations (112) and (115): 
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(118) 


A4 (a = —04) 


£4 47r?7ii 

(118) is independent on the elastic moduli Cij^i which is a consequence of 


^ijk = 0. Thus the dielectric part ^44 of the Green’s function (88) yields 
A4 (a = —04) 1 —1 


^44 (r) = 


(119) 


T 4 + 2 :^ 47r?7ii\/a^p^”+^ 

where 04 = pss/pu- Indeed it is easily checked that (119) is the solution of 
the Poisson equation of a unit point charge (compare (9) and (11)) 


r (V) ^44 (r) + 6^ (r) = 0 
where (compare equation (9)) 

r(V) = 


r (92 

92 ^ 

1 9 "' 


^ 9|/2y 

1 + ''“a? 


Let us consider now the cases / = 1, 2, 3: 

From equations (106)-(110) together with (113) we obtain; 


( 120 ) 


( 121 ) 


Fb (a — —ai) {Cqq — Cii) (C33 — a;C44) + (C13 + C44)^ 

— = A-i = - 


( 122 ) 


Si E, 

These terms coincide with Kroner’s ” A” from [4], Furthermore we obtain: 


Ab± (a — —cii) „ CiiC44af + (Cfg + 2C13C44 — C11C33) ai + C33C44 

- = Dl = - 


Si 


El 


These terms coincide with Kroner’s ”B;” from [4], 
Furthermore we obtain 

Tbc (® ®i) ^ (S44 CIiCqq^ (1-'13 C44) 

— = U; = - 


Si Ez 

These terms coincide with Kroner’s ”C;” from [4], 
Finally we obtain 


(123) 


(124) 
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(125) 


Ac (a — —ai) (C 44 — aiCee) (C 44 — ct,iCn) 

-= - 

These terms concide with Kroner’s ”Dj” from [4], 

Because of the property (116) the elastic part of the Green’s function = 
Gij {i,j = 1,2,3) yields by putting (122)-(125) into (88) and by using the 
abbreviations (89)-(94) Kroner’s result [4]: 


G(r) = E 




X 


( , _ y 2 + ^ 2 ) ^ ^ ^ (^^^2 + 2 ^ 2 ) 

•A-i -;;-h ap A-i -^- 


•Ai 


p* 

xy {aip^ + 2z^] 


y'^z^ — {aip^ + z^) 


xz 

\ 


p^ 


Q- 


Al) 

Cz 


9 ’ 
p2 


P^ 


C, 


xz \ 

p2 

yz 


+ Bp CA— 

p2 


D, 


(126) 


Here we used Kroner’s notation. Above constants B;, C;, D; are defined 
in equations (122)-(125), respectively. 


6 Conclusion 

The electroelastic 4x4 Green’s function of a piezoelectric hexagonal medium 
which is inhnitely extended has been calculated explicitely by using residue 
theory. The obtained expression (88) is highly convenient for a treatment of 
the electroelastic Eshelby tensor [16]. 

An important future application may be the following: The obtained Green’s 
function will be an important quantity even for a treatment of nonlinear, e.g. 
hysteretic material behavior in piezoelectric ceramics. To model these effects 
has become of great interest since the increasing technological importance of 
piezoelectric materials, in particular as electromechanical actuators and sen¬ 
sors [17, 18, 19]. Key contributions for a treatment of such nonlinear effects 
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using the technique of Green’s functions had been presented by Wunderlin 
and Haken [20, 21], 
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